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Abstract: The precise spectra of Cosmic Ray (CR) electrons and positrons have been published by the measurement 
of AMS-02. It is reasonable to regard the difference between the electron and positron spectra ( A<I> = tE>g- —$£+ ) as 
being dominated by primary electrons. The resulting electron spectrum shows no sign of spectral softening above 20 
GeV, which is in contrast with the prediction of the standard model of CR propagation. In this work, we generalize 
the analytic one-dimensional two-halo model of diffusion to a three-dimensional realistic calculation by implementing 
spatial variant diffusion coefficients in the DRAGON package. As a result, we can reproduce the spectral hardening 
of protons observed by several experiments, and predict an excess of high energy primary electrons which agrees with 
the measurement reasonably well. Unlike the break spectrum obtained for protons, the model calculation predicts 
a smooth electron excess and thus slightly over-predicts the flux from tens of GeV to lOOGeV. To understand this 
issue, further experimental and theoretical studies are necessary. 


1 Introduction 

With unprecedentedly high precision, the AMS02 col¬ 
laboration has recently confirmed the excess of positron 
fraction previously observed by PAMELA [J, Q. The 
data shows a discrepancy from the prediction of the stan¬ 
dard model of cosmic ray (CR) propagation (Conven¬ 
tional Model or CM), where the number of events dra¬ 
matically rises at energies ^ 10 — 300 GeV, which is con¬ 
cluded to be an excess of positrons. Proposed explana¬ 
tions include dark matter annihilations/decay ex¬ 

tra sources such as pulsars [1,@], production of secondary 
leptons occuring around acceleration sources [ 3 , 0 |, or 
interactions between high-energy CRs and background 
photons ili^. The precise measurement of AMS also 
gives us an opportunity to study the spectrum of primary 
electrons. Dark matter annihilation/decay and nearby 
pulsar acceleration would produce exccess components of 
positrons and electrons in equal quantities, but because 
the number produced in secondary production by cosmic 
ray interaction with gases is orders of magnitude smaller 
than the number of primary electrons, it is reasonable 
to treat the difference between the flux of electrons and 
positrons as arising from the flux of primary electrons 
only ( Ad) = $6- — d)e+ ). As is pointed out by 11|, 


the spectrum of primary electrons obtained in this way 
exhibits a roughly constant spectral index for energies 
above ^ 20 GeV, while the standard model predicts a 
continuously softening spectrum. This anomaly suggests 
an excess of primary electrons. In order to give a good 
fit to the data, a two-break injection spectrum was found 
inevitable with the allowed parameter space in the CM 

iMl. 

On the other hand, a remarkable spectral harden¬ 


cosmic rays [17 
shock of SNRs 



ing of primary CR nuclei and protons at several hun¬ 
dred GeV has been revealed by ATIC [l^ , CREAM [T^ , 
and PAMELA [T^. Possible explanations of this spec¬ 
tral hardening include the spatial dependent diffusion of 
dual acceleration mechanism at the 
19j| , interaction of cosmic rays with shock 
wave [20| , additional contribution from nearby supernova 
remnants (SNRs) 2l|, the effect of re-acceleration of cos¬ 
mic rays by weak shocks associated with old supernova 
remnants in the Galaxy (23| and so on. 

As assumed in the CM, CRs undergo a uniform dif¬ 
fusion in the galaxy, which provides a good but possibly 
oversimplified approximation of CR propagation. The 
CR diffusion is due to their random scattering by hydro- 
magnetic waves and depends on magnetic field irregular¬ 
ities. On the Galactic disk, the irregularities are caused 
by SNR explosions, while in the diffusive halo, the irreg¬ 
ularities are generated by CRs themselves in the absence 
of SNRs. It has been found that the spectra of the turbu¬ 
lences in the far outer galaxy and in the halo should be 
flatter than that in the local or inner galaxy (2^, which 
implies strong spatial changes of CR diffusion proper¬ 
ties. In this work, we study the impact of latitudinal 
variation of CR diffusion properties on the spectrum of 
primary electrons, and generalize the Two-Halo Model 
(THM) 0 to a three-dimensional realistic calculation 
by implementing spatial variant diffusion coefficients in 
DRAGON package [1^. For comparison, we also per¬ 
form a CM calculation with DRAGON, which assumes 
a uniform diffusion of CRs in space. First, we reproduce 
the spectral hardening of primary protons to determine 
the relevant parameters, and then we attempt to study 
the excess of electrons with this scenario. 
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2 Modeling and parameter settings 


themselves. The diffusion coefficient is expressed as 


It has been long considered that SNRs are the origin 
of primary Galactic Cosmic Rays, and the diffusive shock 
acceleration is regard as the main acceleration mecha¬ 
nism. CRs are accelerated at SNRs and then diffuse 
in the Galaxy, suffering from fragmentation and energy 
losses in the interstellar medium (ISM) and interstellar 
radiation field (ISRF) and magnetic field, decay and pos¬ 
sible reacceleration or convection. Considering those pro¬ 
cesses, the CRs propagation equation can be written as 


dtp{f,p,t) 

dt 


= Q{r,p,t) + '\/ -{D^^Vip-Vctp) 


■0 0 
Tf Tr 


( 1 ) 


where 0(r,p,t) is the density of CR particles per unit 
momentum p at position r, Q{r,p,t) is the source distri¬ 
bution, is the spatial diffusion coefficient, 14 is the 
convection velocity, Dpp is the diffusive reacceleration co¬ 
efficient in momentum space, p = ^ is momentum loss 
rate, and Tf and are the characteristic time scales for 
fragmentation and radioactive decay respectively. In the 
CM, CR diffusion is assumed to be uniform in space and 
only energy-dependent , and the diffusion coefficient is 
parametrized as = l3Do{p/poY, a function where p 
is the rigidity and 5 reflects the property of the ISM 
turbulence. The reacceleration can be described by the 
diffusion in momentum space and the momentum diffu¬ 
sion coefficient is coupled with the spatial diffusion 


coefficient as [251 


D D = _ ^ _ 

pp "" 3d(4-52)(4-d)w 


( 2 ) 


where Va is the Alfven speed, and w is the ratio of mag¬ 
netohydrodynamic wave energy density to magnetic field 
energy density, which can be fixed to 1. In this work, we 
consider a space-dependent diffusion process, and take 
the parameter 6 of the diffusion coefficient as a variation 
along the z axis with d{z). 

To illustrate the effect of spatial change on CR diffu¬ 
sion properties, we first consider a simplified situation. If 
we only consider the diffusion process during CRs prop¬ 
agating in t he g alaxy and neglect other processes, as is 
described in [13 , it can be solved using ID analytical cal¬ 
culations. The proposed THM considers the latitudinal 
changes of CR diffusion properties , and it divides the 
diffusive halo with a typical half-thickness L ~ 5 kpc into 
two regions, where the inner halo is influenced mainly by 
SNRs, with a half-thickness ~ 0.1), and the outer 
halo dominates the rest of the wide region driven by CRs 


D(z for |2;| <^L( inner halo ) 

[Do/ 3(^)"+^, for |z| > ^L{ outer halo ) 

_ (3) 

where p is rigidity and Dq specifies its normalization at 
the reference rigidity Pq. This model behaves like a reser¬ 
voir, where CRs leak out rapidly into the outer halo and 
can re-enter the inner halo. As the result, the spectrum 
of CR primary species is read as 

No = N{z = 0) ~ + (1 -(4) 

which obviously indicates two components of the spec¬ 
trum. We can derive from Eq. o that the inner halo 
dominates the high-energy part with the factor and the 
outer halo mainly contributes CRs at low energies with 
the factor (1—^). Because of the relatively larger diffusion 
coefficient in the outer halo, CRs possess a higher escape 
velocity, which causes fast leakage, and only the low- 
energy CRs re-enter the Galactic disk to the observer’s 
position. Meanwhile, the inner halo with its smaller dif¬ 
fusion coefficient can detain more CRs at high energies. 
As a consequence, we can derive the spectral hardening 
of primary protons from the spatial change of CR diffu¬ 
sion properties. 

In this work, we use the released DRAGON code to 
solve the CR propagation equation described in Eq. ©■ 
DRAGON allows us to perform a numerical calculation 
with a space-dependent coefficient. Without loss of gen¬ 
erality, we take the diffusion coefficient index 6{z) as the 
form 


d{z) 


So + N{-^Y, for l^l <^L( inner halo ) 
(5o-|-A, for l^l >CL{ outer halo ) 


where we use a power-law function to describe the grad¬ 
ual transition from inner halo to outer halo and the index 
n is fixed at 5, representing the extent of how sharply do 
changes to do-|-A. We extend the size of the inner halo to 
with ^ ~ 0.16 to contain the highest number of SNRs. 
For comparison, we carry out calculations for both THM 
and CM using DRAGON, the CM being set with a uni¬ 
form diffusion coefficient in space. Detailed information 
for the parameters of the CR propagation equation is 
shown in Table I. 


Table 1. Some parameters of the propagation equation 


Parameters 

THM 

CM 

Do (cm? s“l) 

5.5x1028 

5.5x1028 

vj^ (kmls) 

25 

25 

Vc (km/s) 

0 

0 

5 

0.19 

0.5 

A 

0.39 
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The injection spectrum of primary CRs at the source 
region is taken as a broken power-law form 


Qi{Ek,r,z) = f,{r,z)qix 


( — ) ““exp (——) 
( — )-“! exp (——) 

'■Pbr' Pcut ' 


if {p<Pbr) 
if {p> Pbr) 
( 6 ) 


where pbr is the break position of rigidity, and p^ut 
represents cutoff rigidity. The normalization condition 
/s(r 0 ,Z 0 ) = 1 is imposed, where Tq = 8.5 kpc is the dis¬ 
tance from the Sun to the Galactic center. The spatial 
distribution of SNRs /s(r,z) is modeled as in [2^. De¬ 
tailed information for the injection spectrum of primary 
protons and electrons is listed in Table 2. 


Table 2. Injection spectrum of primary CRs 


Parameters 

Proton 

Electron 

PbriGV) 

9.0 

5.7 

OiQ 

2.1 

1.7 

ai 

2.31 

2.7 

PcutiGV) 

2.5e6 

20000 


Following the parameters fixed above, we obtain the 
spectrum of primary protons. The result is shown in 

Fig. [II 



Fig. 1. The proton spectra from measurements and 
the two model calculations. The experiment data 
for the proton come from: AMS02 ATICp^, 
PAMELA [3, CREAM [3, KASCADE [13 

The red solid line represents the THM prediction, and 
the black dashed line represents the CM calculation. In 
both calculations, the parameters are tuned to fit well 
with the AMS data from 10 — 10® GeV. At energies 
above ^ 10® GeV, the THM exhibits a pronounced spec¬ 
tral upturn above ~ 10® GeV, and is in good agreement 
with the data within uncertainties. It should be noted 
that although spectral hardening above ~ 10® GeV can 
be recovered, the sharp transition of spectral slope at 
^ 200 GeV found by PAMELA [l^ cannot be repro¬ 
duced. As a matter of fact, when the energy is below 
^ 10 GeV, the solar modulation comes to play a ma¬ 
jor role in the spectrum. The effect of solar modulation 


is taken care of by a so-called force-field approximation 
(^ . with the modulation potential being 550 MV. 

As described above, the effect of spatial change on 
CR diffusion properties is confirmed to produce the spec¬ 
tral hardening of primary protons by a realistic three- 
dimensional model. The observed data above 10® GeV 
support this assumption of a THM. In the next section, 
we will show the result for primary electrons with same 
set of diffusion parameters. 

3 Results of primary electrons 

Unlike primary protons, primary electrons behave in 
a more complicated way. Propagating electrons suffer 
severe energy loss by processes such as synchrotron radi¬ 
ation and inverse Compton scattering, which will rapidly 
steepen the spectrum of electrons with increasing ener¬ 
gies. The power of energy loss is almost proportional 
to A®, so higher energy electrons loss energy rapidly. 
TeV electrons lose most of their energy on a timescale of 
~ 10 ® yr, so only those within a limit of several hundred 
parsecs can reach the solar system. Comparing with the 
power-law or upturn spectra for primary protons in Fig. 
1 , the primary electron spectra are heavily attenuated 
for both THM and CM. The corresponding spectra are 
shown in Fig. 2, together with the subtraction between 
electron and positron flux (A<I> = — $.+ ) measured 

by AMS-02. 



Fig. 2. The electron spectrum in the two models. 
(/)g+ and data are from AMS02 [3^. 13 


In Fig. 2, the observed spectrum seems to be slightly 
upturned at ^ 100 GeV, but there is no sign of a spec¬ 
tral break in the THM (red solid line). The standard 
model (black dashed line) can reproduce the data well 
below ^ 40 GeV, but becomes insufficient as energy in¬ 
creases. The spectrum of primary electrons predicted by 
the THM reveals an integral hardening above ~ 10 GeV, 
which makes it successfully recover the data within error 
bars. However, it can be seen that the THM overesti¬ 
mates the flux of primary electrons slightly with respect 
to the data /\(j) between 30 GeV and 100 GeV. 

Due to rapid energy losses, primary electrons are 
strongly affected by the cooling time, while the influ¬ 
ence of cooling time on primary protons can be neglected 
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above 10 GeV. The cooling timescale of electrons can 
be estimated by Tc ^ 17 Myr (ee/10 GeV)~^ [HI , and 
electrons can diffuse a distance of i? = {2DtcY^^ f^is 
time, where D is the diffusion coefficient. At 10 GeV, 
the distances can be calculated as 3.1 kpc in the GM, 
^2.7 kpc in the inner halo and ^3.2 kpc in the outer halo 
in the THM. It can be seen that the effect of rapid energy 
losses narrows the effective halo size with respect to the 
fixed half-thickness 5 kpc of the halo. The effective halo 
size will continuously shrink with increased energy, so 
the contribution of the inner halo to the total spectrum 
will increase. As described in the former simplified ID 
analysis, the inner halo dominates the high energy com¬ 
ponents with a factor ^ coupled with the contribution of 
the inner halo in the steady state. When the factor ^ 
rises, relatively more high-energy GRs remain in the in¬ 
ner halo and harden the final spectrum. In this scenario, 
the primary electrons predicted by the THM will exhibit 
an integral hardening spectrum above ~ 10 GeV. 

4 Discussion and conclnsion 

In this work, we focus on the latitudinal changes of 
GR diffusion properties and reproduce the spectral hard¬ 
ening of primary protons above ~ 10^ GeV in this sce¬ 
nario. However, the obtained spectrum of primary elec¬ 
trons exhibits an integral hardening with no visible up¬ 


turned break, which gives a good fit with the data above 
~ 100 GeV and overestimates a little between 30 GeV 
and 100 GeV. It should be noted that, although the THM 
produces a non-break spectrum similar to the GM predic¬ 
tion, the additional flux obtained with the THM cannot 
be reached by simply adjusting the global parameters 
used in the GM. The rapid energy losses of GR electrons 
are responsible for the difference between the spectrum 
of primary protons and electrons. 

In addition, it has been pointed out that a continu¬ 
ous source distribution model is not valid for electrons 
in the energy region above 100 GeV [s^, while we still 
consider a continuous distribution in numerical calcula¬ 
tion. The spectrum of electrons at high energies should 
depend on the age and distance of a few nearby sources, 
and some nearby discrete sources have been considered to 
explain the spectrum of electrons 3^, l36| . In the mean¬ 
while, the space-dependent effect of CR diffusion prop¬ 
erties has been used to explain the excess spectra widely 
observed for secondaries l3i1|. th e longitude profile of the 
diffuse 7 — ray emission [l8l |24| , and the significant frac¬ 
tion of the IceCube neutrino flux ]j| . It is worth taking 
both the effect of spatial change on GR propagation and 
nearby discrete sources into account in the future. 

This work is supported by the Natural Sciences Foun¬ 
dation of China (11135010). 
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